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We formulate a stochastic gauge fixing method to study the gauge dependence of Abelian projection. In this 
method, one can change the gauge from the maximally Abelian one to no gauge fixing continuously. We have 
found that the linear part of the heavy quark potential from Abelian contribution depends little on the gauge 
parameter. Similar results have been obtained for the monopole contribution part. We also investigate the gauge 
dependence of the length of monopole loop, which is known to be important for the confinement, and monopole 
density. These results suggest that the picture that monopole plays an important role for the confinement of QCD 
dose not depend on choice of the gauge. 



1. Introduction 

The investigation of quark confinement mecha- 
nism is an important issue in non-perturbative re- 
gion of QCD. As a way to approach the problem, 
Abelian projection has given us many remark- 
able results. For example, Abelian and monopole 
dominances for string tension have been reported 
by G.S.Bali.et.al0, and H.Shiba and T.Suzuki 
have found the infrared effective action of QCD 
in terms of monopole by inverse Monte-Cairo 
method The infrared effective action can 
explain the monopole condensation by energy- 
entropy analysis, and reproduce the string tension 
analytically J3|. 

These results suggest that monopoles which are 
given after Abelian projection have an important 
role for the confinement of the quarks, and sup- 
port the conjecture which has been proposed by 
G.t'Hooft and S.Mandelstam§,|. 

But as is well known, Abelian projection as a 
procedure is explicitly gauge dependent. When 
SU(2) gauge group is reduced into U(l), there 
are infinite ways to fix gauge, and such a remark- 
able results have been reported only in Maximally 
Abelian (MA) gauge. If the conjecture is physi- 
cally meaningful and to understand QCD vacuum 
from the point of view of the monopoles is correct, 
these results should be gauge independent. 
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Recently, Ogilviej6| has developed a charac- 
ter expansion for Abelian and found that gauge 
fixing is unnecessary, i.e., Abelian projection 
yields string tensions of the underlying non- 
Abelian theory even without gauge fixing. Es- 
sentially the same mechanism was observed by 
Ambj0rn and Greensite for Zi center projec- 
tion of SU(2) link variables |?],|| . Furthermore, 
by introducing a gauge fixing function S g f = 
XJ2 ^U fl (x)a3lf fl (x) ji a3, Ogilvie has also shown 
that the Abelian dominance for the string ten- 
sion occurs for small A. Hence he conjectures that 
Abelian dominance is gauge independent and that 
gauge fixing results in producing fat links for Wil- 
son loop and is computationally advantageous for 
the measurements. 

2. Formulation 

At this stage, it is important and significant 
to clarify 1) gauge (in)dependence of monopole 
sector and 2) intermediate region between with 
and without gauge fixing. 

Here we employ stochastic quantization with 
gauge fixing term which has been proposed by 
D.ZwanzigerQ. This method allows us to choose 
gauge fixing condition from exact gauge fixing to 
no gauge fixing continuously. Therefore to apply 
this method to MA gauge fixing yields us a pow- 
crfull tool to research gauge dependence problem 
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of Abelian and monopole dominance. 

For MA gauge fixing, we construct Langevin 
equation for link variables with respect to fictious 
time by analogy with Ref . []lOf : 

U^x^t + St) = w(a;,r) t exp(i/^cr a )C/ M (a;,r) 

XL0(x + fl,T), (1) 
f)C 

uj{x,t) = exp{i——Af at (x,T)cr a 5T). 
2N c a 

In MA gauge, 

Ajot(x,r) = i[a 3 ,X(x,T)] 

= 2(X 2 {x,t)<j 1 -X 1 (x,t)<j 2 ), (2) 

where 

X(x,t) = ^(^(ic, r)(7 3 C^(x s T)t 

A* 

-Up(x - fi,Tya 3 U^(x - /i,t)), 
= ^2Xi(x,T)ai, (3) 

i 

here a — 0(oo) corresponds to exact MA gauge 
(no gauge fixing). 

Then we calculate Wilson loops contributed 
from non-Abelian, Abelian, monopole and pho- 
ton by following way. SU (2) elements can be de- 
composed into diagonal and off-diagonal parts af- 
ter Abelian projection, 



U M (x) = c M (a;)u M (a;), 



(4) 



where c^(x) is the off-diagonal part and u^ix) is 
the diagonal one, 



u^x) = I 
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The diagonal part can be regarded as link vari- 
able of the remaining U(l). One can construct 
monopole currents from field strength of U(l) 
links Q: 

o^(x) = e ll (x) + e v (x + p,)-e M (x + p)-e v (x) 

= Q^v-Vlm^ (-7T < d MV < 7r), 



here, k^(x) is called "monopole current" which 
makes closed loop on 4-dimensional lattice, and 
monopole density is defined by 



X,p, 



(6) 



here iV is a normalization factor. 

Wilson loops from Abelian, monopole and pho- 
ton contributions can be calculated as in Ref. [01 : 
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where d is a lattice forward derivative, d~ is a 
backward derivative and D(x — x') is the lattice 
Coulomb propagator. J„ is the external source 
of electric charge and M^ v has values ±1 on the 
surface inside of Wilson loop. 

We calculate the heavy quark potentials from 
non-Abelian, Abelian, monopole and photon con- 
tributions from these Wilson loops and extract 
the string tensions from each contributions. 

As an improved action to reduce finite lattice 
spacing effects, we adopt the Iwasaki action p2]|. 
We also adopt Runge-Kutta algorithm Jl3|] for im- 
provement of finite step size effects and smearing 
technique for noise reduction. 

3. Numerical Results 

Numerical simulations were performed on 
8 3 x 12, 16 3 x 24 and 24 3 x 32 lattices with 
[3 = 0.995, a = 0.1,0.25,0.5,1.0, and St = 
0.001, 0.005, 0.01. Measurements were done every 
100-1000 Langevin time steps after 5000-50000 
thermalization Langevin time steps. The num- 
bers of Langevin time steps for the thermalization 
were determined by monitoring the functional R 
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Figure 1. Heavy quark potentials from non- 
Abelian, Abelian, monopole and photon contri- 
butions. Lattice size is 16 3 x 24, 5t = 0.005 and 
(3 = 0.995. 



which is maximized in MA gauge and Wilson 
loops. 

In Fig|l] we show the heavy quark potentials 
from Abelian, monopole and photon contribu- 
tions for different a's together with that of non- 
Abelian potential. They can be well fitted by a 
linear and Coulomb terms, respectively. We see 
that the linear parts of potentials are essentially 
same from a = 0.1 to 1.0, and all of them show 
the confinement linear potential behavior. There- 
fore even when we deviate from the MA gauge fix- 
ing condition, we can identify the monopole con- 
tribution of the heavy quark potential showing 
the confinement behavior. As a increases, statis- 
tical error becomes larger. This result suggests 
that the gauge fixing is favorable for decreasing 
numerical errors as pointed by Ogilviejf|. 

In Fig.|| we plot the values of the string ten- 
sions from Abelian, monopole and photon con- 
tributions as a function of the gauge parameter 
a. They are obtained by fitting the data in the 
range 2.0 < R < 7.0. We have taken into account 



string tension 



0.05 



CD 



non-abelian 

{^belian 
I |monopole 
^photon 



Figure 2. String tensions from non-Abelian, 
Abelian, monopole and photon contributions. 
Lattice size is 16 3 x 24, St = 0.005 and (3 = 0.995. 

only statistical errors. The upper two lines stand 
for the range of the non-Abelian string tension. 
The Abelian and the monopole dominances are 
observed for all values of a. On the other hands, 
the string tension from the photon part is con- 
sistent with zero. The string tensions from the 
Abelian parts are about 80% of the non-Abelian 
one. We expect that the difference of the per- 
centage between our result and that of G.S.Bali 
et.al.jjj becomes smaller when we go to larger lat- 
tice size with proper noise reduction technique. 

In FigJ^, we plot monopole densities versus 
gauge parameter a and histograms of monopole 
loops. When gauge fixing condition is far from 
MA gauge, monopole densities has larger value. 
On the microscopic point of view, it has observed 
that short monopole loops which make no cotri- 
bution to string tension increase drastically as a 
increases. Therefore it become difficult to dis- 
tinguish long monopole loop which is responsible 
for reproduction of string tension from short one. 
Then the observables derived from monopole have 
large statistical error at large a region. 

4. Concluding Remarks 

We have developed a stochastic gauge fix- 
ing method which interpolates between the MA 
gauge and no gauge fixing. The method is done 
together with the Iwasaki improved action. 

We have studied the gauge dependence of heavy 
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Figure 3. Monopole densities versus sev- 
eral a(top) and Histograms of short(left) and 
long(right) monopole loops. Lattice size is 24 3 x 
32, St = 0.005 and /3 = 0.995. 

quark potentials derived from Abelian, monopole 
and photon contributions. For Abelian and 
monopole contribution, it is observed that the 
confinement force is essentially independent of the 
gauge parameter. In the calculation of Abelian 
heavy quark potential, we have seen that as gauge 
parameter a increases, the statistical error be- 
comes larger. This result suggests that the gauge 
fixing is favorable for increasing the statistics as 
pointed by OgilvieQ. It is expected that as a 
increase, Abelian string tension would approach 
the non- Abelian one But it is so difficult 

to measure heavy quark potentials in large a re- 
gion, because data are more noisy. In order to 
obtain statistically significant data even in such 
a region of a, more effective noise reduction tech- 
nique such as integral method |L4| will be need. 

For monopole configurations, we have mea- 
sured monopole density and histograms of 
monopole loops. These results show when gauge 
condition is far from exact MA gauge fixing, 
monopole configuration becomes very much com- 
plicated because of increase of short monopole 
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loops and it is difficult to distinguish long 
monopole loop from short one. Therefore in 
large a regions, heavy quark potential from 
monopole contribution has large statistical error. 
If we apply block spin transformation in terms 
of monopole to such a configuration, the clearer 
monoploe configuration would be obtaind. 
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